Geometric phase induced by a cyclically evolving squeezed vacuum reservoir 
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We propose a new way to generate an observable geometric phase by means of a completely 
incoherent phenomenon. We show how to imprint a geometric phase to a system by "adiabatically" 
manipulating the environment with which it interacts. As a specific scheme we analyse a multilevel 
atom interacting with a broad-band squeezed vacuum bosonic bath. As the squeezing parameters are 
smoothly changed in time along a closed loop, the ground state of the system acquires a geometric 
phase. We propose also a scheme to measure such geometric phase by means of a suitable polarization 
detection. 



Whenever a pure quantum state undergoes a parallel 
transport along a closed path, it gathers information on 
the geometric structure of the Hilbert space in which it 
lies. In this letter we will show that a a possible way to 
generate such a parallel transport is by way of an irre- 
versible quantum evolution. In several models of interac- 
tion with the environment there are some "protected" 
subspaces, like the decoherence free subspaces (DFS), 
which are left unaffected Q. States lying in these sub- 
spaces are stationary, i.e. they do not evolve in time. 
A typical example is the ground state of an atomic sys- 
tem, which, trivially, remain unaffected by the coupling 
with the electromagnetic field. However, there are situa- 
tions in which the interaction between a system and an 
engineered environment can generate non-trivial ground 
states 0, S S IE @- F° r instance, when a group of 
atoms collectively interacts with a broad band squeezed 
vacuum, the highly non-classical correlations which are 
present in the field are transferred to the atomic system, 
which relaxes in a complex pure equilibrium state. In 
such a scenario, the control over the engineered reservoir 
allows an indirect control on the state of the system to 
which it is coupled Q ■ Of particular interest is the possi- 
bility to change in time the reservoir parameters in such 
a way that the "protected" system subspace evolves in a 
controlled fashion. Here we show that if this change in 
time is made slowly enough, a state lying in such a sub- 
space evolves coherently and acquires information about 
the geometry of the space explored. 

As an explicit example, we consider a suitable mul- 
tilevel atomic system interacting with a broad band 
squeezed vacuum. To be more specific let us consider 
first a three level atom whose interaction with an elec- 
tromagnetic field in the rotating wave approximation is 
described by the following Hamiltonian: 



H = H s + J tf(u)a(w)du + J [g(u)S f a(u)+H.c.]duj, 
where Hs = SI 2fc=-i k\k){k\ is the free atomic Hamilto- 



nian, S = | — 1) (0| + 10)(1| is the atomic operator describ- 
ing the absorption of an excitation and a(u) is the anni- 
hilation operator of the mode with frequency uj (h = 1). 
The field, which we treat as a reservoir, is assumed to be 
in a broad band squeezed vacuum state. In mathemati- 
cal terms, this is obtained from the ordinary field vacuum 
state by means of the unitary operator K(rj) 



vac (??)) S q = ^0?)|vac), 



(1) 



where 



expi jy [riat (SI - (£1 + w) - H.c] dwj 

(2) 

is a multimode squeezing transformation P, which 
correlates symmetrical pairs of modes around the carrier 
frequency SI and r\ = e %v r is the squeezing parameter, 
whose polar coordinates (p G {0 . . . 2tt} and r > are 
called phase and amplitude of the squeezing, respectively. 

The use of the Born Markov approximation, justified 
by the broadband nature of the field, leads to the fol- 
lowing master equation for the atomic degrees of free- 
dom 0,0: 



dp 



: {R f Rp + pR^R - 2RpR)} 



where T = 2-K\g(Sl)\ 2 , and 

R(t)) = Scoshr + e iv S^ sinhr. 



From (0} follows that the state 
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with c(r) = 



and s(r) 



sinh r 



— . — , — — : — , satisfies 

Vcosn2r v cosh 2r 

R(r])\yjDF(v)) = 0- I n other words, this state is unaf- 
fected by the environment, i.e. it is decoherence free. 
Moreover \ip(r))) represents the new ground state, as all 
the other states of the atomic system relax to it. 
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As anticipated, the key idea is to smoothly change the 
squeezing parameter of the field in order to "adiabati- 
cally" drag a state initially prepared in \tpdf(vo)) into 
\ipDF(Vt)), where r\ t is the time dependent squeezing pa- 
rameter. We will show the existence of an "adiabatic" 
limit such that the transition probability of \iPdf(v)) to 
the orthogonal subspace vanishes as the rate of change 
of rj becomes sufficiently small. Furthermore, we will 
show that after a cyclic evolution of rj, the state \iPdf) 
acquires a geometric phase. It is worth stressing that 
this procedure, although reminiscent of the usual adia- 
batic evolution, is a different physical phenomenon. The 
usual adiabatic approximation refers to a coherent evo- 
lution, obtained by tuning the parameters of the system 
Hamiltonian, while the "steering process" discussed here 
is achieved manipulating the environment. The essential 
difference is that in the latter case the system state can be 
adiabatically controlled entirely by means of an incoher- 
ent phenomenon and no Hamiltonian term contributes 
to its time evolution. To show how this incoherent adi- 
abatic steering process can take place, consider the time 
dependent version of equation J3J where R(r)t) is explic- 
itly dependent on time through rj t ■ It is useful to express 
the equation of motion in the reference frame where tpDF 
is time independent. To this end, consider the following 
unitary transformation 

(c(r)e -1 ^ s(r)e I 2'\ 
1 , (6) 
—s(r)e c{r)e l % / 

from the basis |1), |0), | — 1) to the time dependent basis 
|1), |0), | — 1), where |V>df) coincides with | — 1) Under 
this change of frame, the equation of motion becomes 

^ = -^(&Rp + p&R-2RpR*)-i[G,p\, (7) 

where, in this new frame, p = 0{rj t )pO^ (%), R = 
0{r) t )R(r) t )Oi (r) t ) / Vcosh 2r, f = Tcosh2r and G = 
idO j 'dtO' is a Hamiltonian term arising from the change 
of picture. Moreover, in this frame the Lindbladian term, 
R'R, assumes a simple diagonal form: 

flt J R=|l){l| + |0}(0|. (8) 

The main advantage of this transformation is that it 
allows to formulate clearly the adiabatic condition, since 
the rate of change of the environment parameters are con- 
tained in the operator G. The limit that we are interested 
in is the one in which the dominant contribution in equa- 
tion JJJ comes from the incoherent terms, i.e. \G\ -C V. 

An interesting case is the one in which the squeez- 
ing amplitude is kept constant while its phase is slowly 
changed from to 2n. This adiabatic evolution can be 
easily achieved by tuning, for example, the carrier fre- 
quency, 20, of the squeezed state slightly off resonance 



from the two photon transition | — 1) «-» |1). By intro- 
ducing this detuning, <5, (assuming 8 << Q), the master 
equation obtained has the form of Eq. © and J3J where 
ip is replaced by (pt = (po + St. Hence, a sufficiently small 
value of 5 determines the required adiabatic evolution. 
Under this condition the operator G assumes the form 

. la Q \ 
G = f , (9) 
z \/3 -a) 

where a = cos ^ 2r and j3 = We show that, when 

ip is small enough, the state | — 1) = \4>df) is adiabati- 
cally decoupled from its orthogonal subspace and a cyclic 
evolution in ip results in a geometric phase acquired by 
\iPdf) depending only, in this case, on the amount of 
squeezing r. Note however that, since the steering pro- 
cess is essentially incoherent, any phase information ac- 
quired by a superposition of tpDF and a state belonging 
to the orthogonal subspace is inevitably lost, as the lat- 
ter is subject to decoherence. The only way to retain 
such information is to consider an auxiliary level \a), un- 
affected by the noise, playing the role of a reference state 
for an interferometric measurement. For simplicity as- 
sume that | a) is unaffected by the environment during 
the whole evolution, and, hence, is time independent. As 
a consequence, the action of the unitary transformation 
O on | a) is trivial, and equation (JJJ remains essentially 
unchanged. 

The whole information about the geometrical phase 
and the coherence retained by the system during its evo- 
lution is then recorded in the phase and amplitude of the 
density matrix term p- a = ( — l\p\a), whose evolution 
is described by the following set of coupled differential 
equations 

P-a = -i{a\Gp\l) =i(ap a +- Pp-a) ^ 

P+a = --p +a -i{l\Gp\a) = -i(f + iaip)p +a -if3p- a ^. 

where p+ a = (l\p\a). Assume that initially the excited 
states |1) and |0) of the system are not populated, hence 
p+ a {0) = and the coherence p- a evolves as 

p_„(t) = p_ o (0) • 

7T ^ s [(A- + ia<p) e A +* - (A+ + lap) e x ~ l ] , 



where A± = =F \ y x + iatip — p 1 . 

In the limit r>(/>we obtain for the coherence p- a 

P-a{t) - p- a (0) ((1 - e) e *9**-«f* + ee -<**t-$(l-2e)t) ( 

(10) 

where e ~ ^- . We are interested in a cyclic evo- 
lution, corresponding to T = 2?:. By retaining only the 
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leading terms in ^ the total evolution at time T is given 
by 

p_ a (T) = p_ a (0)e i ^- a ^, (11) 
where we have substituted T with T/a. Finally, going 
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FIG. 1: (a) Schematic representation of the four system con- 
sidered. The energy gap between states | — 1) and |0) and 
between |0) and jl) is Q. The transitions between these level 
are coupled to the modes a(u) of the reservoir. The reference 
state | a) is decoupled from the reservoir, 
(b) Five- level system, transitions 1 «-> and <-> — 1 are cou- 
pled to modes &i(u>) and 1' <-> and «-» — 1' are coupled to 
modes 02 (w) of the reservoir. 

back to the original frame by means of 0'(T)), the cor- 
responding coherence p^ a (t) = (ipDF\p(t)\a) is given by: 



-i7v(l—a) — a/3 -it^ 



(12) 



r 

(13) 

where \rJj(T)) = -fe (e^^l^/WW) + \a)). 

It is clear from this expression, that in the limit 
£ = ip/T <C 1 the dominant contribution to the time 
evolution is just a phase factor e 4 *, with <fi = — tt(1 — a). 
This proves that in the "adiabatic" approximation, the 
system preserves its coherence. In fact, according to 
equation i|12|) , the amplitude damping of p^ a occurs only 
when we take into account the first order contribution in 
£, which shows an exponential decay rate of the order of 
/3 2 ^. This proves that for small ip the system admits an 
adiabatic limit, in which the subspace Hjjp(t) spanned 
by \ipDF{t)) and |a) is adiabatically decoupled from its 
orthogonal subspace H±(t). For this reason, H£> F (t) is 
decoupled from the effects of the decoherence, which only 
affect states lying in its orthogonal subspace. 

Within this approximation, then, a state prepared in 
the space Hdf(0) is adiabatically transported rigidly in- 
side the evolving subspace Hof^)- As a result of this 
adiabatic steering, when the system is brought back to its 



initial configuration, the coherence p^ a acquires a phase 
<p = — ?r(l — a). This phase can be interpreted as the 
geometric phase accumulated by the state |^i?(t)). By 
using the canonical formula for the Berry phase, it easy 
to see that the geometric phase of \ipDF{t)) is given by 



Xg = 



(ip DF \d\ip 



DF) = 



{lpDF\-^\TpDF)dtp 

dtp 



-7r(l — a) = (f>. 



For example a state initially prepared in |V>(0)) = 
^ (IV'-D-f(O)) + |o)), after closing the loop, evolves into 

p(T) = (l-af3\^) \ij(T)))$(T)\+a(3 2 nZ\ip DF )(i> DF \, 



As expected the value of depends only on the squeez- 
ing, and vanishes as the squeezing tends to zero. More- 
over, notice that the phase <$> is purely geometrical, i.e. 
there is no dynamical contribution arising from an ex- 
isting Hamiltonian, since, in absence of any steering pro- 
cess, the states inside Hdf have a trivial dynamics. This 
makes the measurement of this phase a relatively easy 
task. Usual procedures to measure geometric phases 
make use of suitably designed techniques to eliminate 
dynamical phase contributions, such as spin-echo || or 
parallel transport conditions 9]. In this setup, the ge- 
ometric phase is the only contribution to the phase ac- 
cumulated by \ipdf), and hence, it is straightforward to 
measure by a suitable interferometric setup. 

A simple scheme to measure the geometric phase ob- 
tained by such a steering process can be realized with a 
simple variation of our system. Let us consider the five- 
level atomic system shown in picture ^b). It essentially 
consists of two replicas of the three-level system discussed 
above, with the level |0) in common. The important in- 
gredient is that transitions |0) <-* |1) and | — 1) «-> |0) 
arc coupled with modes of the reservoir which are differ- 
ent from those coupled to the transitions |0) |1') and 
I — 1') <-> |0). A simple way to achieve this, is to choose, 
for example, polarisation selective transitions, say, left- 
circularly polarised modes for the former transitions and 
right-circularly polarised for the latter ones. The com- 
plete Hamiltonian of such system is: 



H = H s + / a\(uj)ai(u)duj 

i=l,2 



(14) 



where H s = Q £fe=-i * (|*) (k\ + \k') (fc'|), and Si = 
I - 1)(0| + |0)(1| and S 2 = \- l')(0| + |0)(l'|, anda^w) is 
the annihilation operator of the mode with the energy uj 
and polarization i 6 {1,2}. Assume broadband squeezed 
vacuum states for the set of modes ai(u>) and modes 
ci2(w) with different squeezing parameters 771 = rie lipi 
and 772 = r 2 e l¥ ' 2 : 



|vac(77i, 77 2 )) S(Z = Kx(rn)K 2 (v2)Wa,c), 



(15) 



where Ki{rii) are the analogous of the operator acting 
on the set of modes a,. Under the same assumptions 
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which lead to equation Q we obtain the master equation: 



dp 
~dt 



J2T {R t RiP + pR i Ri ~ 2RipR i } (16) 



where I\ = 2ir\gi(Q)\ 2 , and Ri(r]i) — Sicoshri + 
e llpi Sj sinhr^. This system admits a two-dimensional 
decoherent-free subspace, spanned by states \ipi) and \ip2) 
whose definition is the analogous of state \iPdf) of equa- 
tion (jSJ). We assume again time dependent squeezing 
parameters ry*, and again we examine the time depen- 
dence of the system in a rotating frame, i.e. a frame 
where the state \i/Ji(t)) appear stationary. This leads to 
the following master equation for the five-level system in 
the rotating frame: 



dp 
~dt 



Rip + pR\R 



[Gi,Pl 
(17) 

where Gi = idO/dijiO^rji, 0(t) being the unitary trans- 
formation producing the change of frame. Assume again, 
for simplicity, that the parameters r% and r 2 are kept con- 
stant and that ipx = ip 2 — <P- Under this assumption, the 
master equation can be exactly solved. The solution is 
analogous to the one obtained for system previously ana- 
lyzed. Suppose that the system is initially prepared in a 
coherent superposition of state |V'i( 7 7i)) and |i/> 2 (^2)), for 
example: |^(0)) = ^ (\MVi)) + At a later 

time one has 



s<P (ft ft 



f t, 



with on 



and fa 



cosh In 
eter ip closes a loop, at t 

gained a phase 



T 



(18) 

mn ^ r< . When the param- 

cosn 2ri L 

= 2ir/p, the coherence has 



X = iri (a 2 - ai) = <p 2 



(19) 



which is the difference between the geometric phases 
(pi = tt(1 — acquired by the states \ipi), respectively. 
As in the previous scheme, the visibility is reduced by 
a factor which is linear in the "adiabatic parameters" 
ip/Ti, which guarantees the existence of the adiabatic 
limit. The advantage of this modified scheme is that the 
value of the geometric phases can be readily measured 
from the polarisation of the light emitted when the sys- 
tem relaxes. Infact, if the value of the squeezing param- 
eters ri is suddenly switched to zero, the states \ipi) are 
no longer decoherece free, and decay to a superposition 
of the ground states | — 1) and | — 1'). This dissipation 
process is accompanied by two photon emissions into the 
reservoir. Due to the structure of the interaction Ijl4(l 
with the reservoir, the photon emitted due to the transi- 
tions |1) — > |0) and |1') — > |0), is polarised according to 



the geometric phase accumulated between \ipi) and \ip2)- 
For example, if d\{uj) and 0,2 (u>) are right and left circu- 
larly polarised modes, respectively, the first dissipation 
process produces the linearly polarised photon: 



\ipt) +e l(vi - ¥,2) |^ 2 ) -» \R) +e l ^ 1 - ¥ ' 2) | J L). 



(20) 



The detection of the polarisation of the emitted photon 
makes possible a direct measurement of of the geometric 
phase. 

We have presented a scheme to generate a geomet- 
ric phase via a completely incoherent control procedure. 
This scheme is conceptually different from the usual co- 
herent adiabatic control. The latter is realized through a 
smooth evolution of suitable Hamiltonians, whereas here, 
the adiabatic steering is the effect of an externally con- 
trolled environment. The phase generated is purely geo- 
metrical, and, therefore, experimentally detectable with- 
out resorting to techniques for the elimination of dynam- 
ical contributions. Due to its very nature, this scheme 
is immune from unwanted environmental effects. More- 
over, like any geometric effects, it presents an inherent 
degree of robustness against uncertainties in the control 
parameters. 
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